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ABSTRACT 

We construct a model in which four dimensional chiral fermions arise 
on the boundaries of a five dimensional lattice with free boundary con- 
ditions in the fifth direction. The physical content is similar to Kaplan's 
model of domain wall fermions, yet the present construction has several 
technical advantages. We discuss some aspects of perturbation theory, as 
well as possible applications of the model both for lattice QCD and for 
the on-going attempts to construct a lattice chiral gauge theory. 
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1. Introduction 



The Electro- Weak sector of the Standard Model is a chiral gauge theory. This 
fact motivated numerous attempts to construct a lattice model whose continuum limit 
will be a chiral gauge theory [1-5]. (See ref. [6] for a review). 

The basic stumbling block which was recognized already in Wilson's original 
work [7] is the doubling problem. Following earlier work by Karsten and Smit [8], 
Nielsen and Ninomiya [9] proved that species doubling (and hence a vector-like spec- 
trum) is an unavoidable consequence of any free fermion lattice hamiltonian which 
satisfies some mild assumptions. Thus, attempts of construct chiral lattice gauge 
theories fall into two classes. One approach is to obtains a chiral spectrum at tree 
level by giving up one of the assumptions necessary of the validity of the no-go theo- 
rem[] [1,2]. Typically, the assumption which is violated is that the dispersion relation 
has a continuous first derivative. In the second approach the quadratic hamiltonian 
has a vector-like spectrum, and one tries to eliminate the doublers dynamically by 
introducing some strong interactions which are more effective for the doublers [3-5]. 
However, until now, no lattice gauge model whose continuum limit is consistent with 
the requirements of Lorentz invariance, unitarity etc. has been shown to have a chiral 
spectrum [10-12]. 

Recently, Kaplan [13] proposed that a four dimensional, chiral lattice gauge theory 
may be constructed if one starts with a theory of massive Dirac fermions in five 
dimensions, provided the fermion mass has the shape of a domain wall. Restricted 
to the fifth direction, the free fermionic hamiltonian has a zero mode with definite 
chirality localized on the domain wall. From the point of view of the four dimensional 
domain wall, this zero mode is a chiral fermion. 

For a range of values of the Dirac mass M and the Wilson parameter r, the 
chiral fermion has no doublers if the fifth direction is infinite. But when the fifth 
direction is taken to be finite a doubler appears [13-15]. For example, if we choose 
periodic boundary conditions an anti-domain wall must appear somewhere, and the 
chiral fermion on the anti-domain wall has the opposite chirality. 

As long as one considers the coupling of the fermions to an external gauge field, 
this gauge field can be taken to be five dimensional. But at the dynamical level, 
the introduction of a five dimensional gauge field requires one to deal with a host of 
highly non-trivial issues before the existence of a non-trivial continuum limit of any 
sort can be established. We will not elaborate on this possibility in this paper. 

1 Historically, the work of Drell, Weinstein and Yankielowicz [1] preceded and, to a large extent, 
motivated the work of refs. [8] and [9]. 
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The approach we adopt is to consider the fifth direction as a sophisticated flavour 
space, and not as a full fledged additional space coordinate. The interacting theory 
is defined by coupling the fermions to a four dimensional gauge field [16,17]. This 
means that the link variables satisfy U^(x, s) = U^(x), fi = 1, ... ,4, independently 
of s, that U$(x, s) = 1 and that the gauge field action is any four dimensional lattice 
action which reduces to the standard continuum action in the classical continuum 
limit. Here x^ are the usual four dimensional coordinates and s denotes the fifth 
direction. In the case of a finite lattice we denote the number of sites in the fifth 
direction by N. 

The first question is, are we dealing with a chiral or a vector-like lattice the- 
ory at tree level. The Nielsen-Ninomiya theorem asserts that, assuming (a) complex 
field formulation, (b) hermiticity of the hamiltonian and (c) continuous first deriva- 
tive of the dispersion relation, the massless spectrum of free fermions on a regular 
lattice must consist of an equal number of positive helicity ("right handed") and 
negative helicity ("left handed") Weyl fermion. (One can replace assumption (c) by 
the sufficient conditions that (cl) the number of fermionic degrees of freedom per 
(four dimensional) site is finite and (c2) the hamiltonian has a short range. These 
conditions ensure that the only possible singularities in the dispersion relation are 
level crossings). Moreover, (d) given any charge which is exactly conserved, locally 
defined and has discrete eigenvalues, the equality of the number of left handed and 
right handed fermions holds in every charged sector separately. 

At the level of a free fermion theory, the relevant charges in Kaplan's model 
are the generators of the non-abelian symmetry that we intend to gauge. If the 
fifth direction is infinite, the tree level spectrum can be chiral because assumption 
(cl) is violated. Specifically, for a fixed four momentum one has an effective one 
dimensional lattice problem, whose spectrum consists of a continuous part as well as 
possible bound states. When the energy of the bound state reaches the continuum 
threshold, the bound state disappears from the spectrum, giving rise to a singularity 
in the dispersion relation [13-16]. 

Narayanan and Neuberger [16] considered the introduction of gauge fields di- 
rectly in the "infinite fifth direction setting". To set up perturbation theory, they 
define a four dimensional current J^x) = J2^-oojv( x > s )- Here j^x^s) stands for 
the first four components of the five dimensional current. Apart from an unimportant 
substruction related to the presence of an infinite number of heavy flavours, pertur- 
bation theory is defined in the usual way in terms of correlators of any finite number 
of currents. 

As long as one is interested in perturbative results one can work directly with 
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the Feynmann rules of ref. [16]. But, in order to properly define the model at the 
non-perturbative level, one must consider a sequence of lattice theories with a finite 
fifth direction such that, in the limit N — > oo, the Feynmann rules of ref. [16] are 
reproduced in weak coupling perturbation theory. Since condition (cl) in now fulfilled, 
if the charged fermions spectrum is to remain chiral in these finite lattice models, 
another assumption of the Nielsen- Ninomiya theorem must be violated. Assumptions 
(a) and (c2) are inherent properties of the domain wall model. Consequently, either 
hermiticity of the hamiltonian or tree level gauge invariance (or both) will have to 
be sacrificed. Either way unitarity is violated, as pointed out already by the authors 
of ref. [16]. For this approach to succeed one must construct a specific sequence of 
lattice models with the above properties, and show that all violations of unitarity 
tend to zero in the limit (provided the spectrum is anomaly free in the usual sense). 

If one insists on hermiticity and gauge invariance of the finite lattice action, all 
the conditions of the Nielsen-Ninomiya theorem are fulfilled and so the lattice theory 
is vector-like at tree level [13,14,17]. The question then is whether one can decouple 
the doublers dynamically and achieve a chiral spectrum in the continuum limit of the 
interacting theory. The hope is that the specifics of this model will be sufficiently 
different from previous models which have failed to produce a chiral spectrum [3-5]. 
Whether or not this is the case is the subject of on-going investigations. 

A characteristic property of all versions of the domain wall model is that as the 
number of sites in the fifth direction grows, all correlation functions tend to a limit 
like e~ N . The limiting behaviour is therefore achieved for relatively small N. But it 
remains to be seen whether a consistent chiral gauge theory can be achieved in the 
continuum limit of some version of the model. 

In this paper we have little to say regarding the difficult question of maintaining 
simultaneously a chiral spectrum and a consistent, interacting gauge theory in the 
continuum limit. Our purpose is to eliminate some unessential technical complications 
present in Kaplan's original model, as well as to continue the study of the model at 
the level of perturbation theory. 

In the first part of this paper we construct a five dimensional model in which four 
dimensional chiral fermions arise on the boundaries of a five dimensional slab with 
free boundary conditions in the fifth direction. As in the domain wall model, if the 
fifth direction is (semi)-infinite there is a single chiral fermion on the four dimensional 
boundary. But when the fifth direction is taken to be finite, a doubler appears on the 
other boundary. 

It is well known that the boundary of a sample can have its own dynamics, 
described by an effective field theory in one less dimension. In particular, one finds 
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an intimate relation between Chern-Simons terms in odd dimensions and anomalies in 
even dimensions [18,19]. This relation plays a central role in Kaplan's model [13,20]. 
These ideas have also found interesting applications in the context of the Quantum 
Hall Effect [21]. 

The boundary fermions model is really a variant of Kaplan's model. In fact, a 
careful examination of refs. [13-15] reveals that the zero mode's spectrum is always 
determined by requiring normalizability of the wave function on the side where M and 
r have the same sign. On the other side of the domain wall this condition is always 
fulfilled. If we "discard" that side of the wall we arrive at the boundary fermion 
model. 

The present construction has several technical advantages which merits its dis- 
cussion separately. Analytical expressions, in particular the propagator, take a much 
simpler form in the boundary fermion scheme. In the domain wall model, the prop- 
agator for an infinite fifth direction was calculated in ref. [16]. Here we give the 
propagator both for the semi-infinite and the finite lattice cases. Also, in numeri- 
cal simulations one should obtain the same accuracy in the boundary fermion model 
by taking half as many sites in the fifth direction. This is simply because the fifth 
direction must extend only on one side of the four dimensional boundary. 

In the second part of this paper we discuss the introduction of a four dimensional 
gauge field. We first consider the use of the model for lattice QCD. When slightly 
modified, the model contains an additional parameter m which in QCD plays the role 
of the current mass. Unlike the case of Wilson fermions, we show that perturbative 
correction to the quark mass are proportional to m. Consequently, chiral perturbation 
theory is valid, and the model can be used to study chiral symmetry breaking in lattice 
QCD. 

We next discuss a "mirror fermion" model. In order to construct the model we 
introduce two five dimensional slabs, one for the charged fermions and one for the 
neutral fermions. We also introduce a charged scalar (Higgs) field. We show that in 
the broken symmetry phase, the model can naturally lead to a large hierarchy between 
mirror fermions and ordinary fermions masses. However, there are difficulties, and it 
is not clear whether the model can have a continuum limit describing an interacting 
chiral theory. 

This paper is organized as follows. In sect. 2 we discuss the boundary fermion 
model on a semi-infinite lattice. In sect. 3 we discuss the model on a finite lattice. In 
sect. 4 we consider the introduction of gauge fields. Sect. 5 contains some concluding 
remarks regarding the prospects of obtaining an interacting chiral gauge theory in 
the continuum limit of domain wall or boundary fermions models. 



5 



2. Boundary fermions on a semi-infinite lattice 

We begin with the free field theory defined on a semi-infinite five dimensional 
lattice. For simplicity we work in five dimensions, but everything generalizes to other 
odd dimensions as well. We mainly work with four dimensional momentum eigenstate 
and disregard possible finite size effects in the usual four dimensions. The properties 
of the boundary fermion model presented in this section are similar to those of the 
domain wall model [13-16], and we rederive them here for the convenience of the 
reader. Of course, technical detail such as the explicit form of the bound states and 
the propagator are different (and simpler) in the boundary fermion model. 

the action is given by £ = £4 + £5 where £4 is the usual four dimensional Wilson 
action summed over all s, and £ 5 contains the couplings in the fifth direction. The 
lattice spacing is taken to be a = 1. Explicitly 

f 

£ 4 = E ^( x ' s hi* d id>( x > s) + M E *l>(x, s)i[>(x, s) + -C w , (1) 

x,s,/j, x,s 

£\V = E S ) V M S ) • ( 2 ) 

x,s,/i 

r is the Wilson parameter. We will mainly work with r = 1 and only briefly discuss 
the case r ^ 1. The lattice difference operators are defined by 

df,ip(x, s) = ^{ip(x + fi, s) - i)(x - ji, s)) , (3) 

V M ip(x, s) = ip(x + ft, s) + tp(x — £i,s) — 2ip(x, s) . (4) 
£5 splits into a sum over all s > denoted £' 5 and a boundary term £5, where 

£5 = E s )75<9 5 ^(a:, s) + - E if>(x, s) V 5 ip(x, s) , (5) 

x,s>0 x,s>0 

£5 = 1) + \ 0)GK*» 1) - 21>(x, 0)) . (6) 

^ X ^ X 

Notice that £5 is obtained from eq. (5) by setting s = and dropping terms containing 
fields at the non-existing sites with s = — 1. 

We now go to momentum eigenstates and set r = 1. The action becomes 

£> = I E <M-p> S ) D ^ s '; ?#(p, s ') > (7) 

J P s,s> 

where the integral is over the Brillouin zone, and 

D(s,s';p) = 9(s)0(s')D o (s,s';p). (8) 
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Here 9(s) = 1 for s > and 9(s) = for s < 0. D (s, s';p) is the Dirac operator on 
an infinite fifth direction 

D (s, s'; p) = 1(1 + 7 5 )^+i )S ' + ^ (! ~ 7s)^-i, s ' - (&(p) + , (9) 

p M = sinp M , 

6(p) = l-M + ^(l-cos(^)). (10) 

We now take M to lie in the interval < M < 2. (Later we will see that one 
must further restrict M to lie in the interval < M < 1). The physical content of 
the model is more transparent if we consider the lattice hamiltonian 

H(s,s';p k ) = ^D(s,s';pk,P4 = 0) . (11) 

Notice that in the hamiltonian framework, the fifth direction is really a forth space 
coordinate with a semi-infinite range. We will nevertheless keep calling it the "fifth 
direction". We will show that the spectrum of H(s,s';p k ) contains a right handed 
Weyl fermion that lives on the space boundary. The Weyl fermion is described by 
a bound state with energy E 2 = p 2 ,, which H(s,s';pk) admits for every p k in the 
domain \b(p k )\ < 1. Notice that, for our choice of M, the domain in which chiral 
fermions exist contains the origin of the Brillouin zone but no points in which some 
of the momentum components are equal to ir. As a result, the chiral fermion has no 
doublers [13]. (This situation will change when we make the fifth direction finite). 
The bound state wave function has the form 

* R (s,p k ) = U(s)(l + l5 )^(p k ), (12) 

where the right handed spinor ip(p k ) is a helicity eigenstate 

3 

X! a iPi ^(P*) = E ^(Pk) ■ ( 13 ) 
j'=i 

Substituting this into the eigenvalue equation R = R we find that U (s) must 
satisfy U(s + 1) = b(p k )U(s). The normalized solution is 

U(s) = (l-b 2 (p k ))h s (p k ). (14) 

As promised, this solution is normalizable provided |&(y>fc)| < 1. As |6(pfc)| approaches 
one, the solution decreases slower and slower, until at b{p k ) = 1 it becomes a contin- 
uum eigenstate with vanishing fifth component of momentum. 

Before we turn to the propagator let us briefly discuss the continuous spectrum. 
Because of complete reflection at the boundary, the continuum eigenstates are stand- 
ing waves in the fifth direction. Denoting the fifth component of the momentum by 
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p 5 , the energy of a continuum eigenstate is E 2 = (b(pk) — cosp 5 ) 2 + p\ + p\- If the 
three-momentum pk is fixed, the minimal energy is obtained for P5 = and is given 
by min \E\ = ((b(pk) — l) 2 +Pt) 1 ^ 2 - We see that min \E\ coincides with the energy of 
the bound state on the surface define by b(pk) = 1, which is recognized as the bound- 
ary of the domain in which the chiral fermion exists. The bound state disappears 
when its energy reaches the continuum threshold. This completes our discussion of 
the spectrum of the lattice hamiltonian. 

We now return to the Dirac operator (8) of the eculidean formulation. Unlike 
the Hamiltonian, the eigenstates of the euclidean Dirac operator do not have simple 
physical interpretation. The fact that the Dirac operator is complex allows funny 
things to happen. In particular, the euclidean Dirac operator has a bound state only 
for p M = [22]. On the other hand, the second order operators DD^ and D^D are 
hermitian and non-negative, and so their spectrum is perfectly well behaved. When 
we speak about a bound state spectrum in euclidean space we will always refer to one 
of the second order operators. 

Inspite of the peculiar spectrum of D(s, s';p), the fermionic propagator Gf does 
not show any unexpected behaviour, because Gp — D^G where G is the propagator 
of the second order operator DD^. As a first stage towards the construction of G F , 
we consider the Dirac operator D defined on an infinite s direction (eq. (9)). Going 
to the second order operator DqDq, it is easy to check that the two homogeneous 
solutions are given by exp(±as), where a is defined by the positive solution of the 
equation [16] 

2cosha( P )= 1 + y' +P ' 2 . (15) 

Notice that if b(p) has a zero then a(p) has a logarithmic singularity. To avoid such 
singularities we must take < M < 1. (There is in fact a second allowed range, 
given in In + 1 dimensions by 4n + 1 < M < An + 2. In this case the chiral fermion 
occurs near the corner (rather than near the origin) of the Brillouin zone [15]. The 
physics in both cases is the same, and so we always assume < M < 1). The r.h.s. of 
eq. (15) is always greater that two. This ensures that a(p) is an analytic function in 
the entire Brillouin zone. Notice that both b(p) and e" a(p) tend to 1 - M for p^ -> 0. 
The inverse G of the second order operator -Do-Do is given by 

G (s,s') = Be-^ s - s \ (16) 

B- 1 = 2bsinha. (17) 
The two chiralities are decoupled in the second order operators. In particular 

Drf = \{l + 75)^+ + \{l - 75)H- , (18) 



where Q + and fi_ carry no Dirac indices. This implies a similar decomposition for 
the propagator 

G = i(l + 75)G + + i(l- 7 5)G_. (19) 

We begin with the construction of G_(s, s'). To this end, we apply fi_(s", s) to 
Go(s, s') and check what is the deviation of the result from 5 s » tS >. Since DqDq and f2_ 
contain only nearest neighbour coupling, the deviation vanishes unless s" = 0. For 
s" = one has 

«-(0, s) G (s, s') - V = Be- as '(be- a - 1) . (20) 

s>0 

Similarly, for the homogeneous solutions we find 

5:fi.(^^ = (°; 1 ^JJ' (2i) 

s >o 1 &eT -1, s = 0. 

Eqs. (20) and (21) suggest that G-(s,s') has the form 

G_(s, s') = G (s, s') + K(s') e~ as . (22) 

That is to say, each column of G_ can be constructed from the corresponding col- 
umn of Go, plus a column vector proportional to a homogeneous solution. In order 
that G-(s,s') satisfy physical boundary conditions we allow only the exponentially 
decreasing solution on the r.h.s. of eq. (22). Since the second order operators (and 
hence their propagators) are symmetric, we must have K(s') = A_ exp(-as') where 
the amplitude A_ depends on p^. The amplitude A_ is easily found using eqs. (20) 
and (21). Following the same steps we construct also G + . The final result is 

G±(s,s') =G {s,s') + A ± e- a ^ s+s '\ s,s'>0, (23) 

p a — h 

A- = Be-*>±-±, (24) 

A+ = -Be~ 2a . (25) 

Let us now consider the physical content of these equations. While the amplitudes 
B and A + are regular for all values of p M , the amplitude A_ is singular for small p M , 

„ , . M(2-M) , , x 

■A—ip) = 2 ^ re g u l ar terms , p^ — > . (26) 

This singularity reflects the existence of a bound state whose eigenvalue tends to zero 
for p^ — > 0. The eigenvalue Aq of this bound state can be read off the propagator 
as follows. Using eq. (14) for small p^, the contribution of the bound state to the 
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propagator is approximated by Aq 2 M(2 — M)(l — M)~ s ~ s ' . Equating this to the 
second term on the r.h.s. of eq. (23) and using eq. (26) we find Ag = p 2 as expected. 
Although this time the answer was known beforehand, we have explained this trick 
for extracting a vanishingly small eigenvalue from the propagator because it will be 
useful to us again later. 

We finally notice that the 1/p 2 singularity in G- gives rise to a chiral pole in the 
fermionic propagator 

G F = D ] G 

= -(1 + 7s) M(y2 ~ M \ l- M)- s - s ' + regular terms . (27) 
2 f 

As already mentioned in the introduction, the dispersion relation is singular on the 
surface b{p) = 1 where the bound state disappears from the spectrum. Interestingly, 
the propagator does not show any singularity at b{p) = 1. This is because the 
normalization factor of the bound state tends to zero as b{p) approaches one, and 
so the contribution of the bound state to the propagator vanishes for fixed s and 
s'. One must check, however, whether or not some sort of singularity reappears in 
perturbation theory when summations over the fifth direction are carried out. 

3. Boundary fermions on a finite lattice 

We now proceed to discuss the model with a finite fifth direction < s < N. The 
crucial difference is that now a second chiral fermion appears on the new boundary 
at s = N. Not surprisingly, this fermion has the opposite chirality from the one 
at s = 0. Strictly speaking, there is a tiny mixing between the two chiral modes, 
which vanishes like (1 — M) N . Such exponentially small modifications will cause us 
no concern and henceforth we neglect them. 

The propagator can be constructed using the same method as before. This is 
convenient because the information about the low energy excitations is encoded in 
the singularities of the propagator, and, in any event, what is needed to develop 
perturbation theory is the propagator. 

Denoting quantities that belong to the finite lattice model by a hat, the Dirac 
operator is now given by 

D(s, s') = 9(N - s)6(N - s')D(s, s') . (28) 

Before we actually construct the propagator, it is useful to consider a certain gener- 
alization of the model. We notice that if we start from the Dirac operator defined on 
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a finite fifth direction which has the topology of a circle, then D can be obtained by 
cutting the link connecting the sites s = and s = N. Now, when M is constant, 
using a fifth direction with the topology of a circle (or, equivalently, imposing periodic 
boundary conditions) gives rise to no low energy excitations. Therefore, if starting 
from the Dirac operator of eq. (28) we gradually turn on the link connecting the sites 
s = and s = N, we expect that the two Weyl fermions will form a Dirac fermion 
whose mass is proportional to the strength of this link. 
We are therefore lead to consider the Dirac operator 

in 777, 

D(s, s'; m) = D(s, s') + -(1 - k)S s ,iS s >,n + y (1 + k)S s ,nS s >,i ■ (29) 

As we will see, up to a constant, m indeed plays the role of a Dirac mass for the 
light fermions. Moreover, unlike the case of Wislon fermions, perturbative correction 
will always be proportional to m, i.e. chiral perturbation theory is valid. Thus, the 
model can be used to study chiral symmetry breaking in lattice QCD. (Whether this 
advantage merits the extra trouble involved in going to a five dimensional setting is 
a practical question that we will not address here). 

We now proceed to construct the propagator Gp of the Dirac operator D for a 
general value of m. The propagator for the massless case will be obtained by simply 
setting m = 0. As before, Gf = D^G where G is the propagator of the second order 
operator DD^ . The two chiralities again decouple in the second order operators. 
Furthermore, we now have 

fl + (s,s')=fl-(N-s,N-s'). (30) 

Thus, G+(s, s') is obtained from G-(s, s') by the replacement s,s' N — s,N — s'. 
The notation is the same as in eqs. (18) and (19). 
The propagator must have the following form 

G_(s, s ') = G (s, s') + A_ e - a(s+s "> + A + e -^N-s-s') 

+A m (e- a(N+s - s,) + e -^+*'-*)) . (31) 

This time, the deviation of J2 S DD\s" ', s)G (s, s') from 5 s » iS / vanishes except for s" = 
and s" = N. The same is true for J2 S DD\s", s) exp (±as). Thus, in order to 
construct the s'-th column of the propagator we need a linear combination of Gq(s, s') 
and both of the homogeneous solutions. Taking into account the symmetry of the 
propagator we arrive at eq. (31). 

In solving for the s'-th column of the propagator, the (s'-dependent) coefficients 
of exp (±as) are two unknowns which are determined by solving a two by two matrix 
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equation. There are actually two such equations in which A m appears twice, because 
the s'- dependence of the coefficients can be either exp (as') or exp (—as'). Explicitly 

(32) 
(33) 




where 



the solutions are 



where 



, be° • ur - 1 -mb . 
" 1 -m& 6e a 1 " 



i_ = A- 1 B(l-m 2 )(e a -b), (35) 

i+ = A- 1 fi(l -m 2 )(e- a -6), (36) 

A m = 2A- 1 Bbm cosha , (37) 

A = ft-Met C = e a (6e a - 1) + m 2 (e a - b) . (38) 



In the limit m = we find A m = 0, A_ = A_ and A + = A + . Notice that C 
is diagonal for m — 0, which means that the contributions to the propagator from 
the two boundaries are decoupled. (Strictly speaking, at m = one is left with 
exponentially small off-diagonal terms in C, which can be ignored as long as itself 
is not exponentially small. The infinite (four dimensional) volume limit can always 
be taken in such a way that exponentially small four momenta never occur. All that 
is needed is to take iV to be slightly bigger that the logarithm of the number of sites 
in the ordinary directions). 

The interpretation of m as current mass suggests that we should consider the 
limit where both m 2 and p 2 are small in lattice units. Again, the only amplitude 
which is singular in this limit is We can use the method described earlier to 
extract the smallest eigenvalue from the singular part of the propagator. The result 
is 

|A | 2 =p 2 + m 2 M 2 (2- Mf . (39) 

The current mass of the light Dirac fermion is therefore mM(2 — M). As promised, 
it is proportional to m. 
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4. The interacting theory 



We now proceed to discuss the interacting theory. We first consider the vector-like 
model, which consists of a single five dimensional slab of charged fermions, coupled 
to a four dimensional gauge field as described in the introduction. This model can be 
used to describe lattice QCD, and it poses no conceptual difficulties. Its interesting 
feature is that the current mass m gets only multiplicative renormalization. The 
chiral limit is therefore achieved by letting m tend to zero. In particular, we expect 
that the pion will be massless in this limit. 

In order to verify this picture we have to show that in the limit m — > the Dirac 
fermion remains massless to all orders in weak coupling perturbation theory. This, in 
turn, is true provided the inverse fermion propagator calculated up to n-th order 

rW (s, s'; p) = D(s, s'; p) + gW> (s, s';p) + ...+ g 2n ^ (s, s'; p) , (40) 

has one zero mode on each boundary for p^ = 0. Here £( fe )(s, s';p) is the k-th order 
self energy. 

The physical reason for the stability of the zero mode is the following. In weak 
coupling perturbation theory, corrections to the tree level inverse fermion propagator 
are small. Moreover, since the five dimensional fermion is massive, these correction 
are exponentially suppressed as \s — s'\ grows. Thus, except for an exponentially small 
effect the zero modes living on the two boundaries cannot mix with each other, and 
so a mass term cannot develop. 

A discussion of the perturbative stability of the (single) zero mode has been 
given for the infinite lattice case in ref. [16]. There, the stability was an immediate 
consequence of the absence of additional low energy states that could mix with the 
zero mode. Here we extend the analysis to the finite lattice case and show that, thanks 
to exponential suppression of all correlations in the fifth direction, the masslessness 
of the light fermions is maintained in spite of the fact that the light spectrum is 
vector-like. 

Let us now examine this issue in some detail. The gauge boson propagator is pro- 
portional to S SjS /, and so the one loop self energy is diagonal in s-space SW(s, s';p) = 
S SjS 'T:^(p). The Dirac operator preserves its tree level structure and, for small p M , the 
only change is in the five dimensional mass M — > = M + g 2 YP^{jp = 0). As long 
as < < 1 all the qualitative statements regarding the tree level spectrum and 
tree level propagator apply. In particular, one chiral fermion exists on each boundary. 

At two or higher loop level there appear non-diagonal contributions coming from 
intermediate states of three or more massive fermions. At the n-th order these contri- 
butions decay at least as fast as (1 — M (n ~ 1 )) 3 l s ~ s 'l We have to show that T^(s : s'; 0) 
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has a homogeneous solution which decreases exponentially for s <C N. This is a suffi- 
cient condition for the existence of an (approximate) bound state with exponentially 
small energy near the boundary s = 0. (A similar statement applies for the other 
boundary). 

In order to see that the existence of such a homogeneous solution is stable against 
any small modification of the tree level Dirac operator, it is instructive to first check 
what happens if we take the Wilson parameter r to be close to, yet different from 
1. Consider the Dirac operator defined on an infinite s-direction for r ~ 1. One can 
easily check that for = this Dirac operator has two homogeneous solution for 
each chirality. In the positive chirality sector, one solution behaves approximately 
like (1 — M) s . This is the only solution which is present for r = I. But for r ^ 1 
there is a second solution which behaves like (jEffi) ■ When we restrict our system 
to a semi-infinite s-direction, the boundary term in the Dirac operator picks a linear 
combination of the two solutions. Since both solutions are normalizable for s > 0, 
the result is a normalizable solution too. As r — > 1 the second solution tends to zero, 
until at r = 1 we are left only with the solution (1 — M) s . 

A similar mechanism works for r( n )(s, s'; 0). Imagine that we turn on the elements 
of r( n )(s, s';0) one diagonal at a time. If we consider only s- values which are far 
from both boundaries, the introduction of every new diagonal gives rise to a new 
homogeneous solution. If the new diagonal is below the main diagonal, the new 
homogeneous solution decreases for s > and so it causes no problem. If the new 
diagonal is above the main diagonal, the new homogeneous solution increases for 
s > 0. But when we look for a linear combination of homogeneous solutions which 
satisfies the equation £ s , T (n) (s, s'; O)tf(s') = near s = 0, we find that the boundary 
terms in T^ n \s,s';0) allow us an additional free parameter. With the extra free 
parameter we can find a linear combination which does not include the unwanted 
new solution. 

This completes our discussion of the perturbative stability. Of course, at the non- 
perturbative level we expect that chiral symmetry breaking will take place, giving rise 
to the usual spectrum of confining theories including in particular a massless pion in 
the limit m — > 0. 

We now proceed to discuss the mirror fermion model. This model in obtained by 
taking two five dimensional slabs, each described at tree level by the Dirac operator 
of eq. (29). We next introduce gauge fields only on one slab. This slab will describe 
the charged fermions. The other slab describes neutral fermions. The final step 
involves the introduction of a charged scalar field and the addition to the action of 
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a Yukawa term which couples the boundary layer s = N of the charged fermions to 
the boundary layer s = of the neutral fermions 

C Yukawa = V £ (<j>(x)$(x, N)^°( X , 0) + h.C. ) . (41) 

x 

Here if)(x, s) denotes the charged fermions and ip°(x, s) denotes the neutral fermions. 
The light spectrum consists of one left handed and one right handed charged fermions 
which we denote iPr(x) and as well as one left handed and one right handed 

neutral fermions denoted ip° R {x) and ip° L (x). The right handed fermions arise from the 
two s = boundaries and the left handed fermions arise from the s = N boundaries. 
We repeat this construction for every irreducible representation of the gauge group. 
In the broken symmetry phase we consider the Dirac fermions 




The mass matrix is 



(43) 





Here v is the Higgs VEV. The mass matrix M. is evidently of the seesaw type. Fur- 
thermore, perturbative stability implies that the off-diagonal terms proportional to 
m are only multiplicatively renormalized as in the QCD case. (The effect of quantum 
correction on the Higgs VEV is more subtle. See e.g. ref. [23]). In the limit m <C v, 
the mass matrix M. describes a heavy fermion of mass yv and a light fermion of mass 
m 2 /yv. In particular, taking the limit m — > we find that the massless spectrum 
contains a charged right handed fermion and a neutral left handed fermion. 



5. Prospects 

The massless spectrum of the above mirror fermion model is therefore chiral, 
but the trouble is that in order to achieve this we had to break the gauge symmetry 
spontaneously. If no special fine tuning is made, the gauge bosons mass and the Higgs 
mass will be of the same order of magnitude as the mirror fermions mass. Thus, in 
the low energy limit we obtain a theory of non-interacting chiral fermions. 

The main difference between the model described above and mirror fermion mod- 
els [4] based on Wilson fermions is that no fine tuning is needed in order to obtain 
an undoubled spectrum of light fermions. But the real question is whether one can 
find a version of this model such that at some point in the phase diagram one has 
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simultaneously light gauge bosons and light chiral fermions. Whether models based 
on domain wall or boundary fermions can do better in this regard than models based 
on Wilson fermions is not clear. 

Certain difficulties pertaining to lattice models containing Higgs fields were pointed 
out by Banks and Dabholkar [23]. Some of them can be dealt with by going to a model 
of the Eichten-Preskill type [5] or to a mixed model. Other difficulties are relevant 
even for composite Higgs and are therefore generic. 

Most disturbing is the fact that in all the models discussed in the literature [3-5] 
there does not seem to be a clear mechanism that will distinguish between anomalous 
and non-anomalous theories. This is true even for the Eichten-Preskill model [5]. In 
this model one makes sure that every global symmetry which is broken by instanton 
effects in the continuum model will already be broken explicitly by the lattice action. 
But it is the dynamics which has to determine in what phase the gauge symmetry 
will be realized. 

In the event that the undoubled spectrum is anomalous, the 'tHooft consistency 
condition [24] implies that the theory must be in the Higgs phase. This is because in 
the ungauged model, the non-zero contribution of the light fermions to the anomaly 
can only be cancelled by a Goldstone boson. (The complete spectrum cannot give 
rise to an anomaly because in the underlying lattice theory the charge is exactly 
conserved). But the mechanism which induces the spontaneous breaking of global 
chiral symmetries [25] does not distinguish between would-be anomalous and non- 
anomalous theories. Thus, it is not impossible that spontaneous symmetry breaking 
is a price that must always be paid in order to obtain an undoubled massless spectrum, 
including in theories whose undoubled spectrum is non-anomalous. This concern is 
particularly relevant since in all models one attempts to obtain a chiral massless spec- 
trum before the gauge interactions are turned on. Indeed, recent results [12] provide 
strong evidence that the Eichten-Preskill model undergoes spontaneous symmetry 
breaking, and that the spectrum in vector-like throughout the entire phase diagram. 

The crucial issue is whether in our candidate lattice model there is a mechanism 
which forbids the existence of an interacting chiral continuum limit and which is 
operative only in the absence of anomaly cancellation. In Kaplan's original paper it 
was suggested that such a mechanism does exist in the domain wall model. Namely, 
in the absence of anomaly cancellation there exists a Goldstone- Wilczek current [19] 
away from the domain wall that should prevent the decoupling of the heavy fermionic 
degrees of freedom from the light ones. 

The work of Narayanan and Neuberger [16] represents an attempt to exploit this 
mechanism at the level of perturbation theory. But in order that this mechanism 
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be operative non-perturbatively in a well defined lattice model, we have to sacrifice 
either hermiticity of the hamiltonian or tree level gauge invariance (or both). Typi- 
cally, the resulting unitarity violating effects will arise in a region which is deep indise 
the five dimensional space, far from the region which supports the chiral fermions. 
Consequently, most of the unitarity violating effects will be formally suppressed by 
positive powers of the lattice spacing. (This statement should be true for the effec- 
tive action of the gauge fields obtained by integrating out the fermions). The only 
unsuppressed effect which cannot be cancelled by counter-terms should be the usual 
triangle anomaly. However, even if we arrange for cancellation of the usual triangle 
anomaly, there is a serious danger that additional, finite unitarity violating effects 
will survive and destroy the consistency of the model, because positive powers of the 
lattice spacing can be compensated by divergent loop integrals when we integrate 
over the gauge fields. 

On the other hand, if one insists on hermiticity and gauge invariance of the lattice 
action we do not see how the mechanism described above can be operative. Thus, 
in our view, the prospects of obtaining an interacting chiral gauge theory in the 
continuum limit of models of the kind described above do not seem very promising. 
But more work has to be done before a definite conclusion can be reached. At the very 
least, since these models do not require any fine tuning to maintain the masslessness 
of the light fermions, they can help us focus on the real issue. Namely, the feasibility 
of maintaining simultaneously light gauge bosons and a light chiral spectrum. 
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